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ABSTRACT 

We study the possible contribution of optical galaxies detected with the Hub- 
ble ACS instrument to the near-IR cosmic infrared (CIB) fluctuations in deep 
Spitzer images. The Spitzer data used in this analysis are obtained in the course 
of the GOODS project from which we select four independent 10' x 10' regions 
observed at both 3.6 and 4.5 /im. ACS source catalogs for all of these areas 
are used to construct maps containing only their emissions in the ACS B,V,i,z- 
bands. We find that deep Spitzer data exhibit CIB fluctuations remaining after 
removal of foreground galaxies of a very different clustering pattern at both 3.6 
and 4.5 fim than the ACS galaxies could contribute. We also find that there are 
very good correlations between the ACS galaxies and the removed galaxies in the 
Spitzer maps, but practically no correlations remain with the residual Spitzer 
maps used to identify the CIB fluctuations. These contributions become negli- 
gible on larger scales used to probe the CIB fluctuations arising from clustering. 
This means that the ACS galaxies cannot contribute to the large-scale CIB fluc- 
tuations found in the residual Spitzer data. The absence of their contributions 
also means that the CIB fluctuations arise at z^7.5 as the Lyman break of their 
sources must be redshifted past the longest ACS band, or the fluctuations have 
to originate in the more local but extremely low luminosity galaxies. 
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1. Introduction 

The cosmic infrared background (CIB) is comprised of emission from luminous objects 
from all epochs in the Universe, including by those too faint for individual detection by 
current telescopic studies (see Kashlinsky 2005 for review). A particularly interesting con- 
tributor to CIB comes from the first stars epoch (Santos at al 2002), which may have left 
measurable CIB fluctuations (Kashhnsky et al 2004, Cooray et al 2004). An attempt to mea- 
sure this CIB fluctuations component was made recently by Kashlinsky et al (2005, 2007a,b 
- hereafter KAMMl, KAMM2, KAMM3) using deep IRAC Spitzer images at 3.6 to 8 /xm. 
The first results used a ~ 10 hr exposure of a 5' x 10' field from Fazio et al (2004) and, after 
removing foreground galaxies to uiab ~25 (at 3.6 /xm), detected an excess CIB fiuctuations 
component over that from the remaining "ordinary" galaxies at ^0.5' (KAMMl). KAMMl 
showed that the fiuctuations do not arise from zodiacal light. Galactic cirrus or instrumental 
noise or artifacts (see also KAMM2) and that the residual maps do not correlate with the 
removed sources. A follow-up study by KAMM2 used much deeper GOODS data (Dickin- 
son et al 2003) with ^ 20 hr exposures in two independent sky areas observed at two epochs. 
The new data in four independent fields allowed KAMM2 a better removal of foreground 
galaxies, i.e. until a fixed fioor of the residual shot noise contribution to the power spectrum, 
PsN) was reached. At the same level of Psn all regions have the same large-scale CIB fluc- 
tuations, consistent with their cosmological origin. The new measurements allowed a flrmer 
interpretation of the cosmological sources generating these CIB fluctuations: they must have 
at most a fairly low shot noise component and at the same time produce a signiflcant CIB 
fluctuations component at ^0.5' arising from clustering (KAMM3). 

Recently, Cooray et al (2007) attempted to demonstrate that the signal detected by 
KAMMl at 3.6 /xm originates from optical galaxies detected by the Hubble ACS instrument. 
In one of the GOODS regions they blanked the pixels in the IRAC maps corresponding to 
the optically identifled ACS sources and, from the remaining ~ 20% of the map, evaluated 
the power spectrum of the diffuse emission claiming that it is signiflcantly attenuated when 
the ACS sources are removed. It was pointed out by Kashhnsky (2007) that such analysis 
is flawed: for such deeply cut maps one cannot evaluate clustering properties using Fourier 
transforms because the basis functions are no longer even approximately orthogonal (e.g. 
Gorski 1994). Kashlinsky (2007) further demonstrates that when the correlation function, 
C{6), which is immune to masking, is computed instead for the Cooray et al maps it does 
not depend, within the statistical uncertainties, on whether or not ACS sources are removed. 

In this Letter we present the results from our study of whether the optically detected 
ACS galaxies can explain the CIB fluctuations detected by KAMM. Because the large-scale 
fluctuations do not have a high S/N for the longest IRAC channels, where the instrument 
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noise is significant, we restricted our analysis to the 3.6 and 4.5 fim IRAC channels. We note 
that any explanation of the fluctuations must account for at least the measured amplitude 
and slope at both 3.6 and 4.5 fim and at the right level of the shot noise (KAMM3). We show 
that the galaxies detected in the ACS GOODS survey cannot contribute significantly to the 
GIB fluctuations detected by KAMM because they produce diffuse light with very different 
power spectrum properties and, more importantly, there are negligible correlations between 
the maps used in the KAMM analysis and those of the AGS sources. 



2. IRAC data and results 

The GOODS IRAG data used for this investigation are described in KAMM2. Briefly: 
The GOODS IRAG observations of the HDF-N and GDF-S fields cover regions roughly 
16' X 10' in size. The coverage is obtained through two 10' x 10' mosaics (which overlap 
by ~ 4') that are obtained at epochs ~ 6 months apart (El and E2). For this project the 
individual Basic Galibrated Data (BGD) frames were self-calibrated using the method of 
Fixsen et al (2000), and mosaiced into images with 0'.'6 pixels. The self-calibration solves for 
fixed pattern detector offsets that may be incompletely removed by the BGD pipeline (e.g. 
diffuse stray light), and for frame-to-frame offset variations. The self-calibration does not 
remove the intrinsic large-scale power other than that due to (arbitrary) linear gradients, 
and is particularly useful for studies of the diffuse light fluctuations. Self-calibration cannot 
find a consistent solution if both epochs are combined, because variation in the intensity, 
and especially the direction of the gradient, of the zodiacal light invalidate the assumption 
that the sky itself is a stable calibration source on 6-month timescales. Thus, our analyzed 
regions are four independent sub-fields of 583" x 583" in size with common coverage at both 
3.6 and 4.5 /im. We refer to them below as HDFN-El, HDFN-E2, GDFS-El and GDFS-E2. 

The assembled data were cleaned of resolved sources in two steps: 1) an iterative proce- 
dure was applied whereby each iteration calculates the standard deviation (a) of the image, 
and then masks all pixels exceeding N^nt cr along with A^mask x Ajnask surrounding pixels. The 
procedure is repeated until no pixels exceed Ncut o"- The clipping parameters were fixed so 
that enough pixels ( 1 70%) remain for robust Fourier analysis. As in the early DIRBE work 
(Kashlinsky & Odenwald 2000) the blanked pixels were set to zero, and the resulting power 
spectrum was divided by the fraction of non-zero pixels in the image, thereby preserving the 
total power of the noise. 2) At the second stage we removed a model of the in dividual source s 



described in KAMMl. This is done using a variant of the GLEAN algorithm (iHogbom 19741 ). 
The process identifies the brightest pixel in the image, and then subtracts a scaled PSF at 
that location to remove a fixed fraction of the flux. This is repeated many thousands of 
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times until the components being subtracted are comparable to the noise level of the images. 
Intermediate results are saved after each "iteration" , consisting of the subtraction of ~ 10^ 
components. As progressively fainter sources are removed with each iteration, the shot noise 
from the remaining galaxies decreases as does their contribution to large scale power. 

Comparison of the different regions is most direct if the source cleaning is halted at the 
same specified level of Psn- Fig.l left shows the decrease of Psn with each iteration of the 
cleaning procedure for the four regions used in the study. The right panels in Fig.l show 
the CIB fluctuations derived for the marked level of Psn used in KAMM2. One can see the 
excess over the shot noise at scales ^0.5'. It is important to emphasize that the fluctuations 
have the large scale power spectrum, P{q) oc g^", such that the amplitude of fluctuations is 
approximately flat to slowly rising with scale, corresponding to the effective index n^2. 

3. ACS galetxies and their fluctuations 

We now evaluate the diffuse light fluctuations that ACS galaxies would generate in the 
IRAC maps. For this we constructed synthetic maps over the same regions of the sky at 3.6 
and 4.5 /xm using the HST/ ACS GOODS catalogs (vl.lz; Giavalisco ct al. 2004). These 
were used to generate images of the observed B, V, i, and z-band galaxy distributions that 
are free of (a) small-scale noise and (b) intrinsically large-scale features, whether due to a 
diffuse or unresolved astronomical component or low-level instrumental calibration artifacts. 
The process for each ACS band was as follows: First, four blank images with the same 
scale and orientation of the IRAC images were generated for IIDFN-E1,2 and CDFS-E1,2 
regions. Then, for each cataloged source, the SExtractor (Bertin & Arnouts 1996) parameters 
THETAJMAGE, AJMAGE, and BJMAGE were used to deflne a normalized Gaussian 
shape and orientation for each source, except for "stellar" sources (CLASS_STAR > 0.9) 
which were represented as a single pixel. The intensity of each source was scaled to fit 
the MAG_BEST magnitude, and each source was added to the initially blank image at 
the location designated by the ALPHA_J2000 and DELTA_J2000 parameters. Finally, for 
comparison with the IRAC data, each of these images of the catalog sources was convolved 
with the 3.6 /im or 4.5 /im PSF, as appropriate. For each of the four fields and at both 
3.6 and 4.5 ^m, we prepared one image containing all ACS sources, and four others that 
included only ACS sources fainter than set AB magnitude limits, [mo, mo -|- 2, mo -|- 4, mo -|- 6] 
with mo = [21.8, 21.8, 21.1, 20.6] in the ACS P, V, i, ^-bands. 

After the synthetic maps of ACS sources have been prepared, they were masked using 
the template derived from clipping (Step 1 above) of the IRAC data. Fig. 2 shows the 
fluctuations due to galaxies at all ACS wavelengths in the HDFN-E2 region (others are 
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similar to within statistical uncertainties). Fluctuations of these masked maps show little 
difference out to a certain magnitude, which means that the clipping of IRAC data takes out 
galaxies very accurately for ttiab >23-24 in B, V, i, 2;-band ACS filters. The power spectrum 
of the remaining sources is very different from that of the KAMM CIB fluctuations: it has 
effective index and is very similar to the spectrum of CIB fluctuations detected in 

the deep 2MASS data by Kashlinsky et al (2002), Odenwald et al (2003). The latter arises 
from galaxies fainter than Xvega — 19 (^ab ~ 21 at 2.2 /xm) which are located at 2; ~ 1 
(Cirasuolo et al 2007) as are the bulk of the ACS galaxies. On the other hand the slope of 
the fluctuations in Fig. 1 (right) is described well by the concordance ACDM model at high 
z (KAMM3). The differently sloped power spectra at large scales indicate that the resolved 
ACS galaxies arc unlikely contributors to the large-scale 3.6 and 4.5 //m fluctuations. The 
amplitude of the fluctuations at large scales is also small. Indeed the most that the ACS 
galaxies can contribute to the residual KAMM maps is constrained by the amount of the 
small-scale shot noise power measured to give fluctuations at a few arcsec of order ~ 0.2 
nW/m^/sr at both 3.6 and 4.5 /xm. Scaling the numbers in Fig. 2 by the corresponding 
amount would then give contributions from the remaining ACS sources at arcminute scale 
well below the levels measured by KAMM and shown in Fig. 1 (right) . This by itself implies 
that the KAMM CIB fluctuations cannot arise in the ACS galaxies. 



4. Cross- correlating ACS and IRAC data 

If the ACS galaxies were to explain the KAMM signal there should be a strong correla- 
tion between the ACS source maps and KAMM maps, and their cross-correlation function 
must exhibit the same behavior at scales larger than the IRAC beam. We computed both 
the correlation coefficients, it!o = (<^-^ACS<^-^KAMM)/cACsCKAMM, and the full cross-correlation 
matrix, C{9)^{SFacs{x)SFkamm{x + 0)) between the ACS and KAMM maps. 

Fig 3 (left) shows a very good correlation between the IRAC sources removed by our 
modehng and the sources in the ACS catalog. On the other hand, only very small correlations 
(of order a few %) remain between the residual KAMM maps, containing the fluctuations in 
Fig. 1 (right), and the ACS sources. These correlation coefficients also remain quite small 
as one sub-selects maps made only with progressively fainter ACS sources. 

To further test contributions of ACS galaxies to fluctuations on scales greater than the 
IRAC beam, we computed the correlation function C{9) vs the separation angle 9. In this 
representation, the contributions of any white noise (such as shot noise and/or instrument 
noise) component to C{9) drop off very rapidly outside the beam and for the IRAC 3.6 
lim channel contribute negligibly to the correlation function at 9 greater than a few arcsec. 
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This is best seen from the analog of the correlation coefficient at non-zero lag, defined as 
R{9) = C(^)/(Tacs<^kamm- This quantity is shown in the right panels of Fig. 3 using the 
example of the CDFS-E2 region. The mean square fluctuation on scale 6 is given by the 
integral {6F{6)'^)=-^ C{6')6'd6' (e.g. Landau & Lifshitz 1958). Hence, we evaluated from 
R{6) a related quantity 1Z{6) = ^ R{9')6'd6' , shown in the lower right panels in Fig. 3. 
The correlations arc negligible outside the IRAC beam, which means that, at most, the 
remaining ACS sources contribute to the shot-noise levels in the residual KAMM maps (as 
discussed by KAMM3), but not to the large scale correlations in Fig. 1 (right). 



5. Conclusions 

Our analysis shows that the source-subtracted CIB fluctuations detected recently by 
KAMM in deep Spitzer images cannot originate in the optical galaxies seen in the GOODS 
ACS data. While the B, V, i, z-band galaxies are well correlated with those seen and removed 
in the 3.6 and 4.5 /im data, they correlate poorly with the the residual 3.6 and 4.5 /im 
background emission. These galaxies also exhibit a very different spatial power spectrum 
than the KAMM maps and the amplitude of their fluctuations is generally also low. Thus, as 
discussed in KAMM3, the distant "ordinary" galaxies contribute to the shot noise component 
of the fluctuations, but their contribution to the clustering component is at most small. The 
amplitude of the KAMM source-subtracted CIB fluctuations requires significant CIB fluxes 
from objects fainter than the KAMM subtraction limit. Deep galaxy counts do not show signs 
of turning over at faint magnitudes, but their cumulative CIB levels saturate at magnitudes 
well below the KAMM removal threshold (Madau & Pozzetti 2000, Fazio et al 2004). 

Whatever sources arc responsible for the KAMM fluctuations, they are not present in 
the ACS catalog. There are two ways to reproduce this: 1) Since the ACS galaxies do 
not contribute to the source-subtracted CIB fluctuations, the latter must arise at zt7-5 as 
is required by the Lyman break at rest ~0.1 /xm getting redshifted past the ACS 2;-band 
of peak wavelength ~0.85 //m. This would place the sources producing the KAMM signal 
within the first 0.7 Gyr of the evolution of the Universe and make conclusions of KAMM3 
regarding their very low M/L stronger. 2) Alternatively, the KAMM fluctuations would 
have to originate in lower z galaxies which escaped the ACS GOODS source catalog either 
because they have low surface brightness or are below the catalog flux threshold, but at the 
same time generate diffuse light fluctuations of higher amplitude and different slope than 
the populations already included there. In the latter case we can conservatively estimate 
their expected luminosities, using the z-band mAB~ 24-26.5 corresponding to the source size- 
dependent completeness limit of the ACS catalog (Giavahsco et al 2004). Such galaxies then 
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would have to be very low-luminosity systems since e.g. for mAB=24 the in-band luminosity 
corresponding to {B, z)-band is (4, 2) x lO^/i^^L© at z=l emitted at rest 0.22, 0.45 fim. 

We acknowledge support by NSF AST-0406587 and NASA Spitzer NM0710076 grants. 
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Fig. 1. — KAMM2 results for the 972 x 972 0.6" pixel fields. Triangles correspond to the 
HDFN region, diamonds to CDFS region; open symbols to Ef, filled to E2. Left Show 
decrease of the shot noise contribution to the power spectrum vs the iteration number of 
the KAMM source model. Horizontal lines show the levels reached in KAMM2, which are a 
factor of ~ 2 below those in KAMMl. Right: CIB fluctuations from the residual maps used 
in the present analysis at the shot-noise levels of KAMM2. 
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Fig. 2.— CIB fluctuations due to ACS galaxies for the 972 x 972 0.6" pixel field at HDFN-E2 
region. Filled circles correspond to ACS galaxies fainter than niQ with the mask defined by 
the clipping. Filled diamonds, triangles and squares correspond to fluctuations produced 
by sources fainter than mo + 2, mo + 4, mo + 6. Open diamonds in the 5-band panel show 
the fluctuations produced by galaxies fainter than mAB=23.8 when the clipping mask is not 
applied; the symbols show that such galaxies were effectively removed from the ACS maps 
by clipping alone done in KAMM. Solid hne shows the shot noise slope corresponding to 
the faintest galaxies. It is clear that ACS galaxies are clustered, but the effective index of 
the clustering becomes closer to n ~0 for fainter samples in general agreement with Fig. 2 of 
Kashlinsky et al (2002). 
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Fig. 3. — Left: Correlation coefficient between clipped/masked ACS and KAMM data. Large 
and small symbols correspond to the IRAC Ch 1 and Ch 2. (a) Shows the correlation coef- 
ficient when all ACS catalog sources are included. Open symbols correspond to correlations 
with the maps of the removed sources and filled symbols with the residual KAMM maps 
which contain the fluctuations shown in Fig. 1. Circles, diamonds, triangles, squares cor- 
respond to HDFN-El, HDFN-E2, CDFS-E2, CDFS-E2. In some areas of the sky, such as 
CDFS-E2, a few bright sources (both galaxies and Galactic stars at mAB~21) can notice- 
ably reduce the correlation coefficient, (b) Shows the correlation coefficient with the ACS 
maps which include only sources fainter than the magnitude shown in the horizontal axis. 
Asterisks, pluses, X's and the four-pointed stars correspond to the maps using ACS B,V,i, z 
sources. Right: Black, blue, green, red solid lines show the dimensionless correlation func- 
tion between the diffuse light in the ACS and KAMM maps for B, V, i, z-bands. Dotted line 
shows the dimensionless correlation function of the KAMM maps, C'kamm(6')/o"kamm5 which 
remains positive out to ~ 100" and is better viewed when presented in log-log plots as in 
Fig. SI-4 of KAMMl. The dashed horizontal line marks R = 0. The lower panels show a 
more direct measure of the fluctuations on a given scale, 71(0) = ^ /g R{6')6'd6' . 



